Spin-Orbital Entanglement and Violation of the Goodenough-Kanamori Rules 
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We point out that large composite spin-orbital fluctuations in Mott insulators with t2g orbital 
degeneracy are a manifestation of quantum entanglement of spin and orbital variables. This results 
in a dynamical nature of the spin superexchange interactions, which fluctuate over positive and 
negative values, and leads to an apparent violation of the Goodenough-Kanamori rules. 
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Since the 1960's the magnetism of correlated Mott 
insulators like transition-metal oxides has been under- 
stood by means of the Goodenough-Kanamori (GK) rules 
[ij, ^J. These state that if there is large overlap be- 
tween partly occupied orbitals at two magnetic ions, 
the superexchange interaction between them is strongly 
antiferromagnetic (AF) because of the Pauli principle, 
whereas overlap between partly occupied and unoccupied 
orbitals gives weakly ferromagnetic (FM) interaction due 
to Hund's exchange |2|. In the archetypical case of 180° 
bonds through a single ligand ion this translates into a 
complementary interdependence between spin order and 
orbital order y : ferro orbital (FO) order supports strong 
AF spin order, while alternating orbital (AG) order sup- 
ports weak FM spin order. The canonical example of 
this behavior is KCuFs, where weak FM (positive) spin 
correlations in the ah planes and strong AF (negative) 
correlations along the c axis are accompanied by AG or- 
der in the ab planes and FO order along the c axis. 

The GK rules (and extensions thereof [5j) have been 
extremely successful in explaining the magnetic struc- 
ture in a wide range of materials. This may seem sur- 
prising because they presuppose that the orbital occu- 
pation is static, whereas in recent years it has become 
clear that if partly filled orbitals are degenerate, both 
spin and orbital degrees of freedom should be consid- 
ered as dynamic quantum variables and be described by 
so-called spin-orbital models [y, Q • The GK rules work 
that well because in many compounds a structural phase 
transition, driven by the Jahn- Teller (JT) coupling of de- 
generate orbitals to the lattice, lifts the degeneracy and 
fixes the orbital occupation well above the magnetic tran- 
sition. This happens typically for electrons in Cg orbitals 
where large JT distortions favor C-type orbital order, 
as in KGuFs. However, for t2g orbitals the JT couphng 
is rather weak, and recent experiments in pseudocubic 
perovskite titanates [Sj and vanadates [2| indeed indicate 
that the relevant orbitals fluctuate, and the conditions 
for applying the GK rules are not satisfied. 

In this Letter we investigate the magnetism of corre- 



lated insulators in the case where classical static orbital 
order is absent. We will show that spins and orbitals 
then get entangled due to composite spin-orbital quan- 
tum fluctuations and that the familiar static GK rules are 
violated to the extent that even the signs of the magnetic 
interactions may fluctuate in time. To demonstrate this 
general feature in the most transparent way, we consider 
three different spin-orbital models for correlated insula- 
tors with 180° perovskite bonds between d} , cP and (P 
ionic configurations, respectively, where the GK rules def- 
initely predict complementary signs of spin and orbital 
intersite correlations. The first two models are derived for 
t2g electrons as in LaTiOa (d^) and LaVOs {(f), where we 
demonstrate the violation of the GK rules, while the third 
one is for eg holes as in KCuFa (cP), in which the GK 
rules are perfectly obeyed. This qualitative difference re- 
sults from the quantum nature of t2g orbitals which may 
form singlets, while eg orbitals behave more Ising-like and 
orbital singlets cannot form. 

Superexchange may be regarded to arise from virtual 
excitations into upper Hubbard bands, due to hopping 
with amplitude t, while low-energy charge excitations are 
quenched by strong on-site Coulomb interaction U. The 
resulting spin-orbital models take the generic form 
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where "/ — a,b,c labels the cubic axes. The first term de- 
scribes the superexchange interactions (J — A£^ /U is the 
superexchange constant) between transition metal ions in 
the d" configuration with spin 5*. The orbital operators 
JjY^ and Ki^ depend on Hund's exchange parameter 
rj = Jh /U , which deternfines the spectra of the virtual 
d"d" — > d"^^d"^^ charge excitations. In all three models 
considered here, for each axis 7 only two orbital fiavors 
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are relevant, and J>J' and K^J' can be expressed in terms 

of pseudospin T — 1/2 operators Tj and Tj. Finally, Horh 
stands for the orbital interactions (of strength V) induced 
by the coupling to the lattice — its form depends on the 



type of orbitals {t2g or e^). 

For the <2g systems we will consider chains along the c 
axis, where only two {yz and zx) orbital flavors are active, 
i.e. participate in the hopping. We assume the idealized 
case where these two orbitals contain one electron per 
site, which implies that the third {xy) orbital is empty in 
the d^ model and filled by one electron in the (P model. 
The orbital operators, J^fid^) [IS Ell and J^fid^) [12, 
describing the coupling between the 5' = 1/2 spins of 
the Ti^+ (d^) ions in cubic titanates and that between 
the S = 1 spins of the V"^+ (d^) ions in cubic vanadates, 
respectively, reduce in the absence of Hund's coupling to 
an SU(2)-symmetric expression oc {Ti-Tj + j), which may 
take both positive and negative values. Note also that the 
superexchange [see Eq. (^] thus contains interactions 
like {S+T-){S-T+) + {S-T+){S+T-), which generate 
simultaneous fluctuations of spins and orbitals described 
by the composite operators Qf = SfT~ etc. At finite -q 
both j\f{d^) and Jyf{d'^) also contain 



C^^ = [{{S.rS,)if,-T,))-{S,-S,){f,-T,)]/{2SfA5) 
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This operator appears because double occupancy of ei- 
ther active {yz or zx) orbital is not an eigenstate of the 
on-site Coulomb interaction. Consequently, the total T 
and T^ quantum numbers are not conserved and orbital 
fluctuations are amplified. Finally, GdFeOa-type distor- 
tions induce orbital interactions ex —VTfTj favoring FO 
order along the c axis |l3l| . 

In the Cg system there are two orbital flavors (3z^ — r^ 
and x^ — y^), and for each axis a different linear com- 



bination of them is active (3a: 



along a, 2>y — 
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along 6, and 3z^ — r^ along c). Thus the superexchange 
J-j\d^) between the S = 1/2 spins at the Cu2+ (d^) ions 
in KCuFs is expressed [l^ in terms of axis-dependent or- 
bital operators t[°-'''^ = -\{gI T \/3crf ) and T, 
given by Pauli matrices af and ul- In the absence of 

Hund's coupling 4^V) = (^^^^ " iW]"^ " \). which, 
in sharp contrast to the tig case above, is never negative 
owing to only a single orbital being active along each 
axis. In formal terms, Ji? (d^) is not SU(2)-symmetric, 
and thus orbital singlets are not formed. The Ising-like 
form of J\j{(jP) makes the d^ model look more classi- 
cal than the tig models, but spin-orbital dynamics is still 
promoted as the orbital flavor is not conserved [13. Fi- 
nally, the JT ligand distortions around Cu^+ ions lead to 
orbital interactions ex VT^^ T^"' that favor AO order. 

We investigated intersite spin, orbital and composite 
spin-orbital correlations in the above spin-orbital models. 
To make the results comparable in all cases, we use 



5y = (^ • S,)/{2Sf 
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for the spin correlations. The orbital and spin-orbital 
correlations are defined for the tig {d^ and d^) models as 



while for the Cg (d^) model 
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These definitions of C,V are dictated by the structure 
of the spin-orbital superexchange in the Jh -^ limit. 

We have solved both tig models, d^ and d^, on four- 
site chains along the c axis using periodic boundary con- 
ditions, and we find that nontrivial spin-orbital dynam- 
ics strongly influences the intersite correlations. First 
we consider V^ = 0, i.e. the purely electronic (superex- 
change) spin-orbital models. In the titanate d^ case one 
recovers the SU(4) model 16] in the limit of 77 = 0, 
with robust SU(4) singlet correlations [13 ■ Indeed, in 
the four-site chain all intersite correlations are identical 
and negative, S'y = T^j^ = C^f = -0.25 [Fig. [^a)]. As 
expected, this value is somewhat lower than —0.215 ob- 
tained for the infinite SU(4) chain [l^. At finite 77 one 
finds T-j' < Q- < Sij < as long as the spin singlet 
(5 = 0) ground state persists, i.e. for rj ^ 0.21, and the 
GK rules, which imply that the signs of Sij and T^ are 
different (spin and orbital correlations are complemen- 
tary) are violated. Apparently the composite spin-orbital 

correlations C,v < dominate and cannot be determined 

''J 

from Sij and T^^j by mean-field (MF) decoupling, so the 
spin and orbital variables are entangled, similar to en- 
tanglement in pure spin models J19l |. In fact, the values 
of the correlations indicate that the wavefunction on a 
bond (ij) is close to a singlet of the (total) composite 
quasi-spin Qi + Qj, equivalent to a linear combination 
of (spin-singlet/orbital-triplet) and (spin-triplet/orbital- 
singlet). 

The vanadate d^ model behaves similarly, with all 
three Sij, T^-' and C| ■ correlations being negative in 
the spin-singlet {S — 0) orbital-disordered phase, stable 
for 77 < 0.07 [Fig. db)]. Here the spin correlations are 
weakly AF (5^ ~ -0.05), and AF and FM bonds com- 
pete, promoting a dimerized state [20|. — For both (d^ 
and d^) models the conventional picture is restored at 
large Hund's coupling, which stabilizes the FM ground 
states (at ?7 > 0.21 for d\ and at -q > 0.07 for d^). Here 

the spin-orbital correlations decouple (C|- = 0) and the 
GK rules are perfectly obeyed, with positive Sij = 0.25 
(FM) and negative T^, = —0.5 (AO) correlations. 

The d^ model shows completely different behavior. 
Considering a four-site plaquette in the ab plane, one 
finds that the conventional spin-orbital interrelation 
(AF/FO or FM/AO) is a robust property of the model 
at any value of Hund's coupling. For small 77 ^ 0.25, FO 
correlations T^^ > are accompanied by strong AF spin 



correlations, Sij < 0, and this changes into the oppo- 
site at large rj [Fig. ^c)], just as one would expect from 
the GK rules. Reflecting this situation, the composite 
spin-orbital correlations C^-^ are weaker than the spin 
correlations S^ and the orbital correlations T^^' . 

J tj 



This 



permits spin-orbital separation in the ground state, and 
corrections to this picture are only perturbative |2l|. 

Next we consider finite V, where one expects that the 
coupling to the lattice could suppress the orbital fluctua- 
tions and cure the apparent violation of the GK rules in 
the t2g models. Indeed, at small rj finite V induces orbital 
order and so stabilizes the AF/FO phase [Figs. Cfd) and 
^e)], composite spin-orbital fluctuations are suppressed 
and the GK rules are restored. Already infinitesimal in- 
teraction V > removes the SU(4) symmetry of the (f 
model at 77 = by an Ising-like orbital anisotropy. How- 
ever, for sufficiently large Hund's exchange rj the spin- 
singlet phase survives (unless V :^ J, i.e. orbital in- 
teractions much stronger than the superexchange) . At 
V — J one thus finds three magnetic phases in the d^ {cP) 
model [Figs, dd) and [He)]: (i) AF/FO order ^ in the 
range of 77 ^ 0.04 (77 < 0.06); (ii) an intermediate orbital- 
disordered phase with negative spin, orbital and compos- 
ite spin-orbital correlations of about equal strength, and 
(iii) FM/AO order for 77 > 0.22 (77 > 0.11). The first two 
are separated by an orbital transition within the spin- 
singlet phase. Notably, the GK rules are perfectly obeyed 
in phases (i) and (iii) J23 , while again they do not apply 
in the intermediate phase (ii), which is moved now to a 
more realistic regime of larger 776'. 



In the dp case finite V only stabilizes the large-77 phase 
with FM/AO order at the expense of the small-77 AF/FO 
phase, but the behavior of the model is not changed qual- 
itatively [compare Figs, ^c) and^f)]. We emphasize 
that the different behavior of t2g and eg systems is in- 
trinsic, i.e. has its origin in different spin-orbital physics 
generated by the electronic superexchange interactions. 
In particular, it is not caused by being affected differently 
by coupling to the lattice (i.e., by finite V). 

Further evidence that the GK rules do not directly 
apply in t2g systems follows from the spin exchange con- 
stants Jij = {JiJ ), the expectation value being taken 
over the orbital variables. One finds that in the orbital- 
disordered phase formally FM interaction Jij < is 
in fact, both for ^ = and for finite V, accompa- 
nied by AF spin correlations [Figs. I3a,d) and|2b,e)], 
whereas in the eg case the spin correlations follow the 
sign of Jij for all values of 7] [Fig. I3c,f)]. This re- 
markable difference between t2g and eg systems is due 
to composite spin-orbital fluctuations, which are respon- 
sible for 'dynamical' exchange constants J^J in the for- 
mer case, which exhibit large fluctuations, measured by 
J'ijY^'^ [ll|, as we illustrate here at 
77 = 0. While the average spin exchange constant is 
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J)r fluctuates widely over both positive and 



small in both t 

for d2), J^/) 

negative values. In the d^ case the fluctuations between 

(S = Q/T = 1) and {S = l/T = 0) bond states are so 

large that 5J — 1 \ They survive even quite far from the 

high-symmetry SU(4) point (at rj > 0.1). Also in the d^ 
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FIG. 1: (Color online) Intersite spin Sij (filled circles), or- 
bital T^j (empty circles), and composite spin-orbital Cl^,' 
(crosses) correlations as functions of Hund's exchange rjS, for 
V = (left) and for F = J (right) for: (a,d) d^ model, (b,e) 
d^ model, and (c,f) d^ model. 



FIG. 2: (Color online) Spin exchange constants Jij at y = 
(left) and a.t V = J (right) as functions of Hund's exchange 
775" for (a), (d) d^ model; (b), (e) d^ model; (c), (f) d" model. 
In the shaded regions in (a), (b), (d), (e), Jij is negative (FM) 
and yet the spin correlations are AF, Sij < (see Fig. . 



TABLE I: Energies per site: exact Eq and MF Emf for the 
spin-singlet phases in the three spin-orbital models, obtained 
with four-site clusters. All energies and V in units of J. 
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-0.500 0.0 


-0.316 -0.028 
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-0.655 -0.006 


-0.388 -0.112 
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-0.607 0.097 


-0.311 -0.027 


-0.633 -0.487 



bital correlations that are complementary to each other, 
and both configurations — (orbital-singlet/spin-triplet) 
and (orbital-triplet/spin-singlet) — are entangled in the 
ground state. It remains both an experimental and theo- 
retical challenge to investigate the physical consequences 
of spin-orbital entanglement in real systems. 

This work was supported by the Polish Ministry of 
Science and Education Project No. 1 P03B 068 26. 



model the orbital bond correlations change dynamically 



from singlet to triplet [13| 



resulting in SJ > \Jij\, with 
± 0.247. In contrast, the more 



classical behavior in the dP case is confirmed by 5 J < J, 



as Jij ~ 0.56 and 5. J 



iji 



i{i - (2r,, - \fY/^ c 0.50. 

When quantum entanglement occurs, the ground state 
energy E^ cannot be estimated reliably by MF decou- 
pling of composite correlations (i.e., with the assumption 
C\a ~ 0). The corrections beyond the MF energy Emf 
are largest in the d} case and remain significant in the 
(P model (Table m, but are much less pronounced in the 
<fi model, even ai V — Q. Only when such corrections 
disappear, orbitals disentangle from spins and can be an- 
alyzed separately 24], or spin states can be treated for 
fixed orbital order according to the (static) GK rules. 

We further notice that the (P model exhibits an in- 
teresting property related to the nature of transitions 
between different phases. Namely, the ground state at 
V" = J is a nondegenerate spin-singlet for < 77 ^ 0.11, 
while the orbital quantum numbers change gradually 
from (T) ~ 2 and (T^) ~ ±2 to (T) ~ and (T^) ~ 
in the crossover regime of ry ~ 0.06 [see Fig. ^e)]. We 
have verified that when the orbital terms ex Tl T, are 
neglected, i.e., if Eq. ^ is replaced by an Ising-like term 
T^Tj, a sharp transition occurs instead (from the dou- 
bly degenerate FO state with T^ = ±2 to a disordered 
state with T = 1, T^ = 0), consistent with abrupt tran- 
sitions found before for an infinite chain 25]. Therefore, 
we anticipate that the T^T^ terms induce a continuous 
orbital phase transition in the thermodynamic limit. 

We emphasize that composite spin-orbital fluctuations 
and dynamical exchange constants will control, for realis- 
tic parameters, the behavior of titanates and vanadates. 
In fact, the idea that SU(4)-like fluctuations dominate in 
the ground state has been put forward to understand the 
unusual properties of LaTiOa [ifl and the possible quan- 
tum critical point in the titanate phase diagram [111 |26j . 
Such fluctuations also drive C-AF spin order in LaVOs 
|l2l | and spin-orbital dimerization in YVO3 [Sl, |2a| . 

Summarizing, in correlated insulators with partly filled 
t2g shells, orbitals and spins are entangled, and aver- 
age spin and orbital correlations are typically in conflict 
with the (static) GK rules. These rules should then in- 
stead be understood in terms of dynamical spin and or- 
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